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Abstract 

As a paradigm for heat conduction in 1 dimension, we propose a class of models represented by 
chains of identical cells, each one of which containing an energy storage device called a "tank". 
Energy exchange among tanks is mediated by tracer particles, which are injected at characteristic 
temperatures and rates from heat baths at the two ends of the chain. For stochastic and Hamiltonian 
models of this type, we develop a theory that allows one to derive rigorously - under physically 
natural assumptions - macroscopic equations for quantities related to heat transport, including mean 
energy profiles and tracer densities. Concrete examples are treated for illustration, and the validity 
of the Fourier Law in the present context is discussed. 
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1 Introduction 

Heat conduction in solids has been a subject of intensive study ever since Fourier's pioneering 
work. An interesting issue is the derivation of macroscopic conduction laws from the microscopic 
dynamics describing the solid. A genuinely realistic model of the solid would involve considerations 
of quantum mechanics, radiation and other phenomena. In this paper, we address a simpler set 
of questions, viewing solids that are effectively 1 -dimensional as modeled by chains of classical 
Hamiltonian systems in which heat transport is mediated by tracer particles. Coupling the two 
ends of the chain to unequal tracer-heat reservoirs and allowing the system to settle down to a 
nonequilibrium steady state, we study the distribution of energy, heat flux, and tracer flux in this 
context. 
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We introduce in this paper a class of models that can be seen as an abstraction of certain types 
of mechanical models. These models are simple enough to be amenable to analysis, and complex 
enough to have fairly rich dynamics. They have in common the following basic set of characteristics: 
Each model is made up of an array of identical cells that are linearly ordered. Energy is carried by 
two types of agents: storage receptacles (called "tanks") that are fixed in place, and tracer particles 
that move about. Direct energy exchange is permitted only between tracers and tanks. The two 
ends of the chain are coupled to infinite reservoirs that emit tracer particles at characteristic rates 
and characteristic temperatures; they also absorb those tracers that reach them. To allow for a broad 
range of examples, we do not specify the rules of interaction between tracers and tanks. All the 
rules considered in this paper have a Hamiltonian character, involving the kinetic energy of tracers. 
Formally they may be stochastic or purely dynamical, resulting in what we will refer to as stochastic 
and Hamiltonian models. 

Via the models in this class, we seek to clarify the relation among several aspects of conduction, 
including the role of conservation laws, their relation to the dynamics within individual cells, and 
the notion of "local temperature". We propose a simple recipe for deducing various macroscopic 
profiles from local rules (see Sect. 2.2). Our recipe is generic; it does not depend on specific char- 
acteristics of the system. When the local rules are sufficiently simple, it produces explicit formulas 
that depend on exactly 4 parameters: the temperatures and rates of tracer injection at the left and 
right ends of the chain. 

For demonstration purposes, we carry out this proposed program for a few examples. Our main 
stochastic example, dubbed the "random-halves model", is particularly simple: A clock rings with 
rate proportional to y/x where x is the (kinetic) energy of the tracer; at the clock, energy exchange 
between tracer and tank takes place; and the rule of exchange consists simply of pooling the two 
energies together and randomly dividing - in an unbiased way - the total energy into two parts. 
Our main Hamiltonian example is a variant of the model studied in [14, 9]. Here the role of the 
"tank" is played by a rotating disk nailed down at its center, and stored energy is lo 2 where u is the 
angular velocity of the disk. Explicit formulas for the profiles in question are correctly predicted in 
all examples. 

In terms of methodology, this paper has a theory part and a simulations part. The theory part is 
rigorous in the sense that all points that are not proven are isolated and stated explicitly as "assump- 
tions" (see the next paragraph). It also serves to elucidate the relation between various concepts 
regardless of the extent to which the assumed properties hold. Simulations are used to verify these 
properties for the models considered. 

Our main assumption is in the direction of local thermodynamic equilibrium. For our stochastic 
models, a proof of this property seems within reach though technically involved (see e.g. [4, 19, 8]); 
no known techniques are available for Hamiltonian systems. Extra assumptions pertaining to ergod- 
icity and mixing issues are needed for our Hamiltonian models. It is easy to "improve ergodicity" 
via model design, harder to mathematically eliminate the possibility of all (small) invariant regions. 
In the absence of perfect mixing within cells, actual profiles show small deviations from those pre- 
dicted for the ideal case. 

In summary, we introduce in this paper a relatively tractable class of models that can be seen as 
paradigms for heat conduction, and put forth a program which - under natural assumptions - takes 
one from the microscopic dynamics of a system to its phenomenological laws of conduction. 
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2 Main Ideas 
2.1 General setup 

The models considered in this paper - both stochastic and Hamiltonian - have in common a basic 
set of characteristics which we now describe. 

There is a finite, linearly ordered collection of sites or cells labeled 1, 2, . . . , N. In isolation, i.e., 
when the chain is not in contact with any external heat source, the system is driven by the interaction 
between two distinct types of energy-carrying objects: 

• Objects of the first kind we. fixed in place, and there is exactly one at each site. These objects play 
the role of storage facility, and serve at the same time to mark the energy level at fixed locations. 
For brevity and for lack of a better word, let us call them energy tanks. Each tank holds a finite 
amount of energy at any one point in time; it is not to be confused with an infinite reservoir. We 
will refer to the energy in the tank at site i as the stored energy at site i. 

• The second type of objects are moving particles called tracers. Each tracer carries with it a finite 
amount of energy, and moves from site to site. For definiteness, we assume that from site i, it 
can go only to sites i ± 1. 

With regard to microscopic dynamics, the following is assumed: When a tracer is at site i, it may 
interact - possibly multiple times - with the tank at that site. In each interaction, the two energies 
are pooled together and redistributed, so that energy is conserved in each interaction. The times 
of interaction and manner of redistribution are determined by the microscopic laws of the system, 
which depend solely on conditions within that site. These laws determine also the exit times of the 
tracers and their next locations. A priori there is no limit to how many tracers are allowed at each 
site. We stress that this tracer-tank interaction is the only type of interaction permitted: the tanks 
at different sites can communicate with each other only via the tracers, and the tracers do not "see" 
each other directly. 

All stochastic models considered in this paper are Markovian. Typically in stochastic rules 
of interaction, energy exchanges occur when exponential clocks ring, and energy is redistributed 
according to probability distributions. In Hamiltonian models, tracers are usually embodied by 
real-life moving particles, and energy exchanges usually involve some types of collisions. 

The two ends of the chain above are coupled to two heat baths, which are infinite reservoirs 
emitting tracers at characteristic temperature (and also absorbing them). It is sometimes convenient 
to think of them as located at sites and N + 1. The two baths inject tracers into the system 
according to certain rules (to be described). Tracers at site 1 or N can exit the system; when they 
do so, they are absorbed by the baths. The actions of the two baths are assumed to be independent 
of each other and independent of the state of the chain. The left bath is set at temperature T L ; the 
energies of the tracers it injects into the system are iid with a law depending on the model. These 
tracers are injected at exponential rates, with mean g L . Similarly the bath on the right is set at 
temperature T R and injects tracers into the system at rate g R . 

To allow for a broad spectrum of possibilities, we have deliberately left unspecified (i) the rules 
of interaction between tracers and tanks, and (ii) the coupling to heat baths, i.e., the energy distri- 
bution of the injected tracers. (Readers who wish to see concrete examples immediately can skip 
ahead with no difficulty to Sects. 3.1 and 4.1, where two examples are presented.) We stress that 
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once (i) and (ii) are chosen, and the 4 parameters T L , T R , g L and g R are set, then all is determined: 
the system will evolve on its own, and there is to be no other intervention of any kind. 

Remark 2.1 Our approach can be viewed as that of a grand-canonical ensemble, since we fix the 
rates at which tracers are injected into the system (which indirectly determine the density and energy 
flux at steady state). An alternate setup would be one in which the density of tracers is given, with 
particles being replaced upon exit. In this alternate setup, the 4 natural extensive variables would be 
the temperatures T L and T R , the density of tracers (mean number of tracers per cell) and the mean 
energy flux. For defmiteness, we will adhere to our original formulation. 

We now introduce the quantities of interest. For fixed N, let /j, N denote the invariant measure 
corresponding to the unique steady state of the ./V-chain (assuming there is a unique steady state). 
The word "mean" below refers to averages with respect to fi N . The main quantities of interest in 
this paper are 

• s i = mean stored energy at site i ; 

• e i = mean energy of individual tracers at site % ; 

• k i = mean number of tracers at site i ; 

• E i = mean total energy at site i, including stored energy and the energies of all tracers present. 

For simplicity, we will refer to e i as tracer energy and E i as total-cell energy. 

We are primarily interested in the profiles of these quantities, i.e., in the functions i s^, e i: k i 
and E i as N — > oo with the temperatures and injection rates of the baths held fixed. More 
precisely, we fix T L , T R , g L and g R . Then spacing {1,2,..., N} evenly along the unit interval 
[0,1] and letting N — ► oo, the finite-volume profiles i s i ,e i ,k i ,E i give rise to functions 
£ i—*- s(0, e(£), &;(£), E(£), £ € [0, 1]. It is these functions that we seek to predict given the mi- 
croscopic rules that define a system. 

2.2 Proposed program: from local rules to global profiles 

We fix N, T L , T R , g L and g R , and consider an iV-chain with these parameters. To determine the 
profiles in Sect. 2.1, we distinguish between the following two kinds of information: 

(a) cell-to-cell traffic, and 

(b) statistical information pertaining to the dynamics within individual cells. 

In (a), we regard the cells as black boxes, and observe only what goes in and what comes out. 
Standard arguments balancing energy and tracer fluxes give easily the mean energy and number of 
tracers transported from site to site. While these numbers are indicative of the internal states of the 
cells (for example, high-energy tracers emerging from a cell suggests higher temperatures inside), 
the profiles we seek depend on more intricate relations than these numbers alone would tell us. 

We turn therefore to (b). Our very naive idea is to study a single cell, and to bring to bear on 
chains of arbitrary length the information so obtained. We propose the following plan of action: 

(i) Consider a single cell plugged to two heat baths (one on its left, the other on its right), both 
of which are at temperature T and have injection rate g, T and g being arbitrary. Finding 
the invariant measure /j, T ' e describing the state of the cell in this equilibrium situation is, in 
general, relatively simple compared to finding n N . 
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(ii) Suppose the measure fi T,e has been found. We then look at an iV-chain with T L = T R = T 
and g L = g R = g, and verify that the marginals at site i of the invariant measure fi N are equal 
to fi T,B . (By the marginal at site i, we refer to the measure obtained by integrating out all 
variables pertaining to all sites / i.) 

(iii) Once the family {fi T ' Q } is found and (ii) verified, we assume that the structure common to the 
fi T ' e passes to all marginals of fi N as N — ► oo even when (T L , g L ) / (T R , g R ). More precisely, 
for all £ 6 (0, 1), we assume that all limit points of the marginals of fi N at site [£iV] (where 
[x] denotes the integer part of x) inherit, as N — ► oo, the structure common to fi T,e . 

We observe that (i)-(iii) alone are inadequate for determining the sought-after profiles, for they 
give no information on which T and g are relevant at any given site. The main point of this program 
is that (a) and (b) together is sufficient for uniquely determining the profiles in question. 

Remark 2.2 Our rationale for (iii) is as follows: Fix an integer, i. As N — > oo, the gradients 
of temperature and injection rate on the £ sites centered at [£N] tend to 0, so that the subsystem 
consisting of these £ sites resembles more and more the situation in (ii). Though rather natural from 
the point of view of physics [3], this argument does not constitute a proof. Indeed our program is in 
the direction of proving the existence of well defined Gibbs measures and then assuming, when the 
system is taken out of equilibrium, that thermodynamic equilibrium is attained locally; in particular, 
local temperatures are well defined. The full force of local thermal equilibrium is not needed for 
our purposes; however. The assumption in (iii) pertains only to marginals at single sites. 

The rest of this paper is devoted to illustrating the program outlined above in concrete examples. 

3 Stochastic Models 

3.1 The "random- halves" model 

This is perhaps the simplest stochastic model of the general type described in Sect. 2.1. The mi- 
croscopic laws that govern the dynamics in each cell are as follows: Let 5 > be a fixed number. 
Each tracer is equipped with two independent exponential clocks. Clock 1, which signals the times 
of energy exchanges with the tanks, rings at rate | \fx where x is the (current) energy of the tracer. 
Clock 2, which signals the times of site-to-site movements, rings at rate y/x. The stored energy at 
site i is denoted by y { . In the description below, we assume the tracer is at site i. 

(i) When Clock 1 rings, the energy carried by the tracer and the stored energy at site i are pooled 
together and split randomly. That is to say, the tracer gets p(x + y^ units of energy and the 
tank gets (1 — p)(x + y ri ), where p £ [0, 1] is uniformly distributed and independent of all other 
random variables. 

(ii) When Clock 2 rings, the tracer leaves site i. It jumps with equal probability to sites i ± 1. If 
i = 1 or N, going to sites or N + 1 means the tracer exits the system. 

It remains to specify the coupling to the heat baths. Here it is natural to assume that the energies of 
the emitted tracers are exponentially distributed with means T L and T R . 

This completes the formal description of the model. 

Remark 3.1 The rates of the two clocks are to be understood as follows: We assume the energy 
carried by the tracer is purely kinetic, so that its speed is ^fx. We assume also that a tracer travels, on 
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average, a distance 5 between successive interactions with the tank, and a distance 1 before exiting 
each site. 

Remark 3.2 As we will show, the invariant measure does not depend on the value of 5, which can 
be large or small. The size of 5 does affect the rate of convergence to equilibrium, however. 

Remark 3.3 While the tracers do not "see" each other in the sense that there is no direct interaction, 
their evolutions cannot be decoupled. The number of tracers present at a site varies with time. When 
two or more tracers are present, they interact with the tank whenever their clocks go off, thereby 
sharing information about their energies. A new tracer may enter at some random moment, bringing 
its energy to the pool; just as randomly, a tracer leaves, taking with it the energy it happens to be 
carrying at that time. 

3.2 Single-cell analysis 

3.2.1 Single cell in equilibrium with 2 identical heat baths 

We consider first the following special case of the model described in Sect. 3.1: N = 1, T L = T R = 
T, and Q L = Q R = Q. Each state of the cell in this model is represented by a point in 



where Q, k = . . . , x k }, y) : x#,y £ [0, oo)}. Here {xi, . . . x k } is an unordered /c-tuple 

representing the energies of the k tracers, y denotes the stored energy, and a point in Q k represents 
a state of the cell when exactly k tracers are present. 

Remark 3.4 We motivate our choice of Q. During a time interval when there are exactly k tracers 
in the cell - with no tracers entering or exiting - it makes little difference whether we think of the 
tracers as named, and represent the state of the cell by a point in [0, oo) fc+1 , or if we think of them 
as indistinguishable, and represent the state by a point in Q k . With tracers entering and exiting, 
however, thinking of tracers as named will require that all exiting tracers return later, otherwise the 
system is transient and has no invariant measure. Since any rule that assigns to each departing tracer 
a new tracer to carry its name is necessarily artificial, and for present purposes exact identities of 
tracers play no role, we have opted to regard the tracers as indistinguishable. 

We clarify the relationship between [0, oo) fc+1 and Q k and set some notation: Let ir k : [0, oo) fc+1 
— > Q k be the map 7r k (x±, . . . , x k , y) = ({x 1 , . . . x k }, y), i.e., ir k is the (fc!)-to-l map that forgets the 
order in the ordered fc-tuple (x 1: . . . , x k ). For a measure jl on [0, oo) fe+1 that is symmetric with 
respect to the x £ coordinates, if /x = (vr^/i, and a and a are the densities of jj, and p, respectively, 
then a and a are related by 



We also write d{x l , . . . , x k }dy = (Tr k )^(dx 1 . . . dx k dy), and use / to denote the characteristic 
function. 



oc 



n 




(disjoint union) 



o-({x 1 ,...,x k },y) = kla(x 1 ,...,x k ,y) . 



Proposition 3.5 The model in Sect. 3.1 with N = 1, T L = T R = T, and Ql = = Q has a unique 
invariant probability measure fi = fi T,e on £1 This measure has the following properties: 
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• the number of tracers present is a Poisson random variable with mean x = 2q^Jit /T, i.e., 




fc = 0,l,2,... ; 



(3-1) 



• the conditional density of ' \i on Q k is c k a k d{x 1 , 



,x k }dy where 



a k ({x 1 ,...,x k },y) = I{ Xl ,..., Xk ,y>o} 



1 



e -/3(x 1 +-+x fc +j/) 



y/X 1 ■ ... ■ X k 



(3.2) 



here (3 = 1/T, and c k = (3 k\ (f3/ir) k / 2 is the normalizing constant. 



Proof: Uniqueness is straightforward, since one can go from a neighborhood of any point in Q to a 
neighborhood of any other point via positive measure sets of sample paths. We focus on checking 
the invariance of \x as defined above. 

For z, z' € Q, let P h (dz'\z) denote the transition probabilities for time h > starting from z. 
We fix a small cube A C £l k for some k, and seek to prove that 



On the time interval (0, h), the following three types of events may occur: 

Event E 1 : Entrance of a new tracer 

Event E 2 : Exit of a tracer from the cell 

Event E 3 : Exchange of energy between a tracer and the tank 

We claim that with initial distribution fi, the probability of more than one of these events occurring 
before time h is o(h) as h — ► 0. This assertion applies to events both of the same type and of 
distinct types. It follows primarily from the fact that the these events are independent and occur at 
exponential rates. Of relevance also are the exponential tails of a k and the Poisson distribution of 
p k := /u(J7 fc ) in the definition of fi. To illustrate the arguments involved, we will verify at the end of 
the proof that the probability of two or more tracers exiting on the time interval (0, h) is o(h), but 
let us accept the above assertion for now and go on with the main argument. 

Starting from the initial distribution ji, we let P(E i ,A) denote the probability that E i occurs 
before time h resulting in a state in A, and let V{E{ n E\ n E$, A) denote the probability of starting 
from a state in A and having none of the E i occur before time h. We will prove 



P(£ 1; A) + P(E 2 , A) + P(£ 3 , A) + P(E1 n E c 2 n E c 3 , A) - ^{A) = o{h)^{A) . (3.3) 



Notice that A can be represented as the union of disjoint sets UjAj where each A { is of the form 

A e(z) = {({^D • • -, x k},y) ■ x e G [x e ,x e + £],£= 1, . . . , k,y e [y,y + e]} 

for some z = ({x 1 , . . . , x k },y) G £1, e <C 1, and with the intervals [x e , x t + e] pairwise disjoint 
for t = 1, 2, . . . , k. To prove (3.3) for A, it suffices to prove it for each A { provided o(h) in (3.3) is 
uniformly small for all i. 




h=0 



= 0. 
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We consider from here on A = A e (z) with the properties above. Let a denote the density of jjl. 
With e sufficiently small, we have /i(A) « a(z)e k+1 = p k c k a k (z)e h+l where c k and a k are as in 
the proposition. The other terms in (3.3) are estimated as follows: 

P(E 1 , A): E 1 is in fact the union of 2k subevents, corresponding to a new tracer coming from the 
left or right bath and the k approximate values of energy of the new tracer. For definiteness, we 
assume the new tracer arrives from the left bath, and has energy in [x 1 ,x 1 + e]. That is to say, the 
initial state of the cell is described by 

B = {({x 2 ,...,x k },y) : x e G [x e ,x e + e],y G [y,y + e]} C %_i . 

The contribution to P(E 1 , A) of this subevent is 

fi(B)hg [ Xl+£ Pe~ 0x dx « hg ji{B) e~ m ^ (3e 

J x 1 

~ hg p k _ x c k-\ a k (z) ^JT X (ie k+1 . 

Here, hg is the probability that a tracer is injected, and the integral above is the probability that the 
injected tracer lies in the specified range. Summing over all 2k subevents, we obtain 

k 

P(E 1 , A) « 2hg(3 (J2 y/fy) Pk-ic-k-i e k+1 . (3.4) 
t=\ 

F(E 2 ,A): In order to result in a state in A, the initial state must be in 

c i = {({xi,---,x k ,x},y) : x e G [x e ,x e + e],x G [0,oo),y G [y,y + e]} C U k+1 . 
We assume here that the tracer with energy x exits between time and time h. This gives 

f°° 1 
P(E 2 ,A) « c s+1 a k (z) e k+l / min{hv^, 1} -=e _/J:e dx 

Jo v x 

= P k+1 c k+1 a k (z) e k+1 (5-\h+o{h)) . (3.5) 

P(E 3 , A): For definiteness, we assume it is the tracer with energy near x 1 that is the product of the 
interaction with the tank. To arrive in a state in A £ , one must start from 

D = {({ x i, ■ ■ ■ ,^fe}>y) : x t 6 [x e ,x e + e] for £ > 2,x l +y G [% + +y + 2e]} . 
A simple integration using the rule of interaction in Sect. 3.1 gives 



J>(E 3 ,A) « h^± p- k c k a- k (z)e k+1 . (3.6) 

P(E% f] E 2 f] E 3 , A): We first note that starting from A, the probability of the tracer with energy 
w x x exiting is 



Pk c k <T k^)^/x~ l e Px ^e k / hy/x-=e-P x dx = h^/x\ p k c k a k (z) 

Jx 1 V x 



£ k+l 
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the probability of the tracer with energy « x 1 interacting with the tank is 



and the probability of a new tracer entering the cell from the left (resp. right) bath is hg fi(A £ ). Thus 
PCEf n£ 2 c fl E%, A) « PfcC E a s (z) e S+1 • II^l - /ly 7 ^) • (1 - hg) 2 ■ U e (l - hy/¥ t /S) 

« p m (z) £ 5+1 - (l-/ l (sf =lV ^ + 2 e + isf =lV ^)). (3.7) 

Summing Eqs. (3.4)-(3.7), we obtain (3.3) provided 

2g(3p h _ 1 c ]! _ l = p- k c- k and 2gp s c s = Pk +l c~ k+l T . 
Note that these two equations represent the same relation for different fe. We write this relation as 

c kPk = 

Ck+iPk+i %Q ' (3 ' 8) 
and verify that it is compatible with assertion (3.1): Since 

and J °° x^^e^^dx = VttT, we have 

To complete the proof, we estimate the probability of two or more tracers exiting before time h. 
For n = 2, 3, . . . , let E 2n be the event that the initial state is in 

c n = {({ x i,---,Xk,x (1 \...,x {n) },y) : x e G [x e ,x e + e],x (n G [0,oo),y G [y,y + e]} , 

and during the time interval (0, h), all n of the tracers carrying energies x (£ \ £' = 1, 2, . . . , n, exit 
the system. The probability of U n >2^2,n * s 

Y,Pk + nC- k+ n^)e k+1 (0(h))\ 

n>2 

which, from Eq. (3.8), is bounded by 

Pm(z)e k+1 (^Y(W) n = Pm(z)e k+1 o(h) . 

n>2 ^ ' 



□ 



Remark 3.6 In the setting of Proposition 3.5, since the cell is in equilibrium with the two heat 
baths, it is obvious that it ejects, on average, 2g tracers per unit time, and the energies of the tracers 
ejected have mean T. We observe that the cell in fact reciprocates the action of the bath more 
strongly than this: the distribution of the energies of the ejected tracers is also exponential. To see 
this, fix k and consider one tracer at a time. The probability of the tracer exiting with energy > u is 



r r 1 

Ju V x 



e-^dx = I3~ l e-P u 
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3.2.2 Chain of N cells in equilibrium with 2 identical heat baths 

We treat next the case of arbitrary N. That is to say, the system is as defined in Sect. 3.1, but with 
T L = T R = T and g L = g R = g. Let fi be as in Proposition 3.5. 

Proposition 3.7 The N-fold product fix ■ ■ ■ x fiis invariant. 

Remark 3.8 That the invariant measure is a product tells us that at steady state, there are no spatial 
correlations. We do not find this to be entirely obvious on the intuitive level: one might think 
that above-average energy levels on the left half of the chain may cause the right half to be below 
average; that is evidently not the case. This result should not be confused with the absence of 
space-time correlations. 

Proof: Proceeding as before, we consider a small time interval (0, h), and treat separately the 
individual events that may occur during this period. One of the new events (not relevant in the case 
of a single cell) is the jumping of a tracer from site % ± 1 to site i. We fix a phase point 

z = (#>,..., i = ({xf x«},^>; {4 N \...,4 N M N) )> 

and let A = l\f =l A {i) C Q N where A (i) = A £ {z {i) ) is as in Proposition 3.5. Let a (i) = a k .(z (i) ). For 
definiteness, we fix also an integer n, 1 < n < N, and assume that at time 0, the state of the chain 
is as follows: 

(i) at site n + 1, there are k n+l + 1 tracers the energies of which lie in disjoint intervals 

[4 n+1 \ 4 n+1) + £],..., {x[ n+l \ 4 n+1) + e] and z< n) + e] , 

(ii) at site n, there are k n — 1 tracers whose energies lie in 

[^2 ^ ' -*4 ^ + £ ] ' • • • ' [^4j* ' **4j* £ ] • 
urring and tl 

into site n is then given by \l ■ E ■ M where 



— (fl) 

The probability of this event occurring and the tracer with energy x\ jumping from site n + 1 



J = IX 



^n,n+l (p k C k .O®e k i +1 ) 



fxf+e i 

M = Pk +1 c k +1 a (n+1) e k "+i +1 / h^—e'^dx 



This product can be written as 

h (TIN „ „ «)Jc, + l\ ^n-l^n-l Pfe n + 1 + 1 C *„+l+l 



f (n~ l!W ^«) - — — - — — ■ y , 

K n K n K n+1 K n+1 

which, by Eq. (3.8), is equal to 

\ i^l,P k c k /^).^. (3.9) 
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There are many terms of this kind that contribute to J I A (z')P h (dz'\z)n(dz), two for xf for 
each pair (n, £). We claim that the system has detailed balance, i.e., the term associated with the 
scenario above is balanced by the probability of starting from a state in A and having the tracer at 
site n carrying energy jump to site n + 1. The probability of the latter is 

ldlfL lPk c k <r^ +1 ) ■ Jx^e^- ■ H * h^c^=e~^dx , 
2 ki k > v 1 e J S M y/x 

which balances exactly (3.9) as claimed. 

An argument combining the one above with that in Proposition 3.5 regarding the injection of 
new tracers holds at sites 1 and N. □ 

Propositions 3.5 and 3.7 are steps (i) and (ii) in the proposed scheme in Sect. 2.2. 
3.3 Derivation of equations of macroscopic profiles 

Having found a candidate family of equilibrium measures {p, T,B }, we now complete the rest of 
the program outlined in Sect. 2.2. The next step, according to this program, is to assume that for 
N S> 1, the marginals of the invariant measure p, N at site i are approximately equal to fi T,e for 
some T = Tj and g = g^ We identify those parts of our proposed theory that are not proved in this 
paper and state them precisely as "Assumptions". 

Assumption 1. Given T L ,T R > 0, g L , g R > 0, and N £ Z + , the N -chain defined in 
Sect. 3. 1 with these parameters has an invariant probability measure p, N . 

We do not believe this existence result is hard to prove but prefer not to depart from the main line 
of reasoning in this paper. Once existence is established, uniqueness (or ergodicity) follows easily 
since any invariant measure clearly has strictly positive density everywhere. A proof of the statement 
in Assumption 2 below is more challenging. For £ 6 (0, 1), let Hn^n] denote the marginal of p, N 
at the site [£N]. 

Assumption 2. For every £ G (0, 1), every limit point as N — > oo of A*jv,[£iV] * J a 
member of the family {p 7 ^' 8 , T > 0, g > 0}. 

In Sect. 2.1, we introduced four quantities of interest. There is one that was somewhat ambigu- 
ously defined, namely e t . Its precise meaning is as follows: e i := J^k^iPi,^^ where p i k is the 
probability that the number of tracers at site i is equal to k and e i k is | of the mean total tracer 
energy conditioned on the number of tracers present being equal to k. 

Theorem 3.9 is about the profiles of certain observables as TV —> oo. We refer the reader to the 
end of Sect. 2.1 for the precise meaning of the word "profile" in the theorem. 

Theorem 3.9 The following hold for the "random-halves" model defined in Sect. 3.1 with arbitrary 
T L ,T R , g L , g R . Under Assumptions 1 and 2 above: 

• the profile for the mean number of jumps out of a site per unit time is 



3(0 = 2(g L + (g R -g L )0 
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• the profile for the mean total energy transported out of a site per unit time is 

Q(0 = 2(g L T L + (e R T R -e L T L )0; 

• the profile for the mean stored energy at a site is 



s(0 = 



M) Ql + - 



in the case g L = g R , this simplifies to s(£) = T L + (T R — T L )£ ; 

• the profile for mean tracer energy is e(£) = ^s(0 ; 

• the profile for mean number of tracers is 



"(0 = \I^M)\ 



s(0 



the profile for mean total-cell energy is 



E(0 = s(0 + x(Oe(0 = s(0+ 2^(0 3(0- 
Proof of Theorem 3.9 : We divide the proof into the following three steps: 

/. Information on single cells: Items (i)-(iv) are strictly in the domain of internal cell dynamics. 
The setting is that of Proposition 3.5, and the results below are deduced (in straightforward compu- 
tations) from the invariant density given by that proposition. The parameters are, as usual, T and g. 
(Note that this means the rate at which tracers enter the site is 2g.) 

(i) stored energy has density (5e~^ y and mean T ; 

(ii) tracer energy has density -^J=e~ l3x and mean ; 

(iii) mean number of tracers, x = 2- s Jl^g ; 

(iv) mean total-cell energy, E = T ■ (1 + ^) . 

Items (v) and (vi) are in preparation for the analysis of cell-to-cell traffic : 

(v) mean number of jumps out of the cell per unit time, j = 2g ; 

(vi) mean total energy transported out of the cell per unit time, Q = Tj. 

II. Global phenomenological equations: Consider now a chain with N cells with settings T L ,T R , g L 
and g R at the two ends. The following results use only standard conservation laws together with the 
local rule that when a tracer exits a cell, it has equal probability of going left and right. 

(A) Balancing tracer fluxes: Let j\ denote the number of jumps per unit time out of site i. Then 



^ = 2 [0l + 1 ^ [ (Qr-QO) ■ (3-10) 

Proof: Consider an (imaginary) partition between site i and site i + We let — Aj i denote the flux 
across this partition. Then Ajj = ^(jj+i — 3d for i = 0, 1, . . . , where j and jjv+i are defined 
to be 2g L and 2g R respectively. For i ^ 0, N, the \ is there because only half of the tracers out of 
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site i + 1 jump left, and half of those out of site i jump right. The fluxes across partitions between 
all consecutive sites must be equal, or there would be a pile-up of tracers somewhere. This together 
w i tn J2i A?i = 2(g R — Qi) gives the asserted formula. 

(B) Balancing energy fluxes: Let Q { denote the mean total energy transported out of site % per unit 
time. Then 

Proof: The argument is identical to that in (A), with AQ { = \{Q i+1 — Q { ) and Q and Q^ + i 
defined to be 2g L T L and 2g R T R respectively. 

III. Combining the results from I and II: Fix £ G (0, 1). Passing to a subsequence if necessary, we 
have, by Assumption 2, Mat,kat] fi T ^ ,e ^ for some T(£) and £>(£). We identify the two numbers 
T(0 and q(0 as follows: 

Let j N {0 denote the mean number of jumps out of site [£N] per unit time in the iV-chain. Then 
by (A) in part II, := li m Ar^oo j N (0 = + (Qr ~ This is, therefore, the mean number 

of jumps out of a cell with invariant measure fj T ^ ),e( ^\ Similarly, we deduce that the mean total 
energy transported out of a cell with the same invariant measure is Q(£) := lim^^^ Q N (Q = 
2(q l T l + (q r T r — QiT L )£). Appealing now to the information in part I, we deduce from (v) and 
(vi) that q(0 = i j(0 and T(£) = Q(0/j(0- 

The rest of the profiles follow readily: we read off s = T from (i), and deduce the relation 
between s and e by comparing (i) and (ii). The profiles for x and E follow from (iii) and (iv) 
together with our knowledge of q and T. The proof of Theorem 3.9 is complete. □ 

Remark 3.10 It is instructive to see what Theorem 3.9 says in the special case when g L = 0. Since 
no particles are injected at the left end, clearly, T L cannot matter. But since tracers do exit from 
the left, one expects an energy flux across the system. Upon substituting g L = into the formulas 
above, one gets 

j(0 = 2Q R £ , e(0 = i S (0 = \T R , x(£) = 2g R i^jT R ' , 
and an energy flux of — \Q'{Q = —QrT r . 

3.4 Simulations 

Numerical simulations are used to validate Assumptions 1 and 2 in Theorem 3.9 . 

We mention here only those details of the simulations that differ from the theoretical study. 
Needless to say, we work with a finite number of sites, usually 20. Simulations start in a random 
initial state, and are first run for a period of time to let the system reach its steady state. All times are 
in number of events (E x , E 2 , and E 3 ). Up to half the simulation time is used to reach stationarity; 
statistics are then gathered during the remaining simulation time. Since total simulation time is 
finite, we find it necessary to take measures to deal with tracers of exceptionally low energy: these 
tracers appear to remain in a cell indefinitely, skewing the statistics on the number of tracers. We 
opted to terminate events involving a single tracer at a single site after about 0.001 of total simulation 
time. This was done about 10 times in the course of 10 9 events. 

Simulations are performed both to verify directly the properties of the marginals at individual 
sites and to plot empirically the various profiles of interest. Excellent agreement with predicted 
values is observed in all runs. A sample of the results is shown in Fig. 1. 
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Figure 1: Random-halves model with 20 sites, temperatures T L = 10, T R = 100 and injection rates 
g L = 10, g R = 5. Top left: Mean tank energies s^. Top right: Mean number of tracers >c i . Bottom: 
Mean total energy E { . Simulations in perfect agreement with predictions from theory. 
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3.5 Interpretation of results 

We gather here our main observations, discuss their physical implications, and provide explanations 
for the reasons behind the derived formulas. The content of this subsection is valid not only for the 
random-halves model but for all of the examples studied in this paper. 

1. Linearity of profiles 

We distinguish between the following 3 types of profiles: 

a. Transport of energy and tracers: and Q(Q are always linear by simple conservation laws. 

b. Mean stored and (individual) tracer energies: s(0 and e(£) are linear if and only if there is no 
tracer flux across the system. (See item 4 below.) 

c. Tracer densities and total-cell energy: x(£) is never linear (unless T L = T R and g L = g R ); in 
addition to the obvious bias brought about by different injection rates, tracers have a tendency to 
accumulate at the cold end (see item 3 for elaboration). As a result, E(£) is also never linear. 

2. Heat flux and the Fourier Law 

Heat flux from left to right is given by $ = T L g L — T R g R . Thinking of the temperature of the 
system as given by T(£), Theorem 3.9 says that thermal conductivity is constant and proportional 
to T R — T L if and only if there is no tracer flux across the system, i.e., if and only if g L = g R . 

3. Distribution of tracers along the chain 

In the case g L = g R , more tracers are congregated at the cold end than at the hot. This is because 
the only way to balance the tracer equation is to have the number of jumps out of a site be constant 
along the chain. Inside the cells, however, tracers move more slowly at the cold end, hence they 
jump less frequently, and the only way to maintain the required number of jumps is to have more 
tracers. When g L / g R , the idea above continues to be valid, except that one needs also to take into 
consideration the bias in favor of more tracers at the end where the injection rate is higher. 

4. Tracer flux and concavity of stored energy 

One of the interesting facts that have emerged is that s(£) is linear if and only if g L = g R , 
and when g L / g R , their relative strength is reflected in the concavity of s(£). This may be a 
little perplexing at first because no mechanism is built into the microscopic rules for the tanks to 
recognize the directions of travel of the tracers with which they come into contact. The reason 
behind this phenomenon is, in fact, quite simple: If there is a tracer flux across the system, say from 
right to left, then the tank at site i hears from site i + 1 more frequently than it hears from site i — 1 
(because > It therefore has a greater tendency to equilibrate with the energy level on 

the right than on the left, causing s i to be > ^(s i+1 + Sj_i). Since this happens at every site, a 
curvature for the profile of s { is created. The reader should further note that tracer flux and heat flux 
go in opposite directions if (g L — g R ) ■ (g L T L — g R T R ) < 0. 

5. Individual cells mimicking heat baths 

The cells in our models are clearly not infinite heat reservoirs, yet for large N, they acquire 
some of the characteristics of the heat baths with which they are in contact. More precisely, the ith 
cell injects each of its two neighbors with g i tracers per unit time. These tracers, which have mean 
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energy T ir are distributed according to a law of the same type as that with which tracers are emitted 
from the baths (exponential in the case of the random halves model); see Remark 3.6. Unlike the 
conditions at the two ends, however, the numbers Q i and T i are self-selected. 

3.6 A second example 

We consider here a model of the same type as the "random-halves" model but with different micro- 
scopic rules. The purpose of this exercise is to highlight the role played by these rules and to make 
transparent which part of our scheme is generic. 

The rules for energy exchange in this model simulate a Hamiltonian model in which both tracers 
and tanks have one degree of freedom. Write x = v 2 and y = uj 2 , v,u> € (— oo, oo), and think of 
energy as uniformly distributed on the circle {(v,u>) : v 2 + lo 2 = c}, so that when the tracer interacts 
with stored energy, the redistribution is such that a random point on this circle is chosen with weight 
\v\ (this is the measure induced on a cross-section by the invariant measure of the flow). That is to 
say, if (x, y) are the stored energy and tracer energy before an interaction, and (x', y') afterwards, 
then for a£[0,i + y], 



P{y> > a} = P{\u'\ > M = 1 - J0 V = 1 _ (3 U) 



or, equivalently, the density of y' is 

1 1 



for y' £ [0, x + y] 



Assume now that all is as in Sect. 3.1, except that when Clock 1 of a tracer rings, it exchanges 
energy with the tank according to the rule in (3.11) and not the random-halves rule. Following the 
computation in Sect. 3.2 (details of which are left to the reader), we see that Propositions 3.5 and 
3.7 hold for the present model provided a k is replaced by 



c k ({x 1 , - . . .-,x k },y) — I{ Xl ,...,x k ,y>0} 



-/3(x 1 +...+x fc +{/) 



V x i • • • • • x kV 

This defines a new family of {fi ' 8 } for this model. With fi T ' 8 in hand, we make the assumption 
as before that for N » 1, the marginals of individual sites have the same form. Proceeding as in 
Sect. 3.3, we read off the following information on single cells: 

(i) stored energy has density const. e"^/ ^Jy and mean T/2 ; 

(ii) tracer energy has density const. e~ f3x / \fx and mean T/2 ; 

(iii) mean number of tracers, k = ; 



(iv) mean total-cell energy, E = ^(1 + x) . 

The rest of the analysis, including (iv), (v), (A) and (B), do not depend on the local rules (aside 
from the fact that tracers exiting a cell have equal chance of going left and right). Thus they remain 
unchanged. Reasoning as in Theorem 3.9, we obtain the following: 



Proposition 3.11 Under Assumptions 1 and 2, the profiles for the model with energy exchange rule 
(3.11) are 
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• Q(0 = 2(q l T l + ( % T R - ^T L )0 ; 

• 8(0 = e(0 = \ Q(0/M) ; 

. x(0 = y/2ir/s(OXO/2; 

. E(0 = s (0 + x(Oe(0 = s(0 + V^Ui(0/2. 

Numerical simulations give results in excellent agreement with these theoretical predictions. 
4 Hamiltonian Models 

In Sect. 4.1, we introduce a family of Hamiltonian models generalizing those studied numerically 
in [14, 9]. A single-cell analysis similar to that in Section 3 is carried out for this family in Sect. 4.2, 
and predictions of energy and tracer density profiles are made in Sect. 4.3. We again use the As- 
sumptions in Section 3, but the predictions here are made under an additional ergodicity assumption, 
ergodicity being a property that is easy to arrange in stochastic models but not in Hamiltonian ones. 
Results of simulations for a specific model are shown in Sect. 4.5. A brief discussion of related 
models is given in Sect. 4.6. 

4.1 Rotating disks models 

We describe in this subsection a family of models quite close to those studied numerically in [14, 9]. 
The rules of interaction (though not the coupling to heat baths) are, in fact, used earlier in [17]. 

4.1.1 Dynamics in a closed cell 

We treat first the dynamics within individual cells assuming the cell or box is sealed, i.e., it is not 
connected to its neighbors or to external heat sources. 

Let T C M 2 be a bounded domain with piecewise C 3 boundary. In the interior of T lies a 
(circular) disk D, which we think of as nailed down at its center. This disk rotates freely, carrying 
with it a finite amount of kinetic energy derived from its angular velocity; it will play the role of 
the "energy tank" in Sect. 2.1. The system below describes the free motion of k point particles (i.e., 
tracers) in V = T \ D. When a tracer runs into dT , the boundary of T , the reflection is specular. 
When it hits the rotating disk, the energy exchange is according to the rules introduced in [17, 14, 9]. 
A more precise description of the system follows. 

The phase space of this dynamical system is 

& k = (T k xdDx R 2k+1 )/ ~ 

where 

x = (xi, . . . , x k ) G T fe denotes the positions of the k tracers, 

•& S dD denotes the angular position of a (marked) point on the boundary of the turning disk, 
v = (v l , . . . , v k ) 6 M. 2k denotes the velocities of the k tracers, 
uj G R denotes the angular velocity of the turning disk, 

and ~ is a relation that identifies pairs of points in the collision manifold M k = {(x, v, u) : x t G 
dF for some £}. The rule of identification is given below. 
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The flow on Cl k is denoted by «3> s . As long as no collisions are involved, we have 

$ s (x, v, to) = (x + s\, i? + soj, v, lo) . 

We assume at most one tracer collides with dT at any one point in time. (I> s is not defined at multiple 
collisions, which occur on a set of measure zero.) At the point of impact, i.e., when x t G dT for one 
of the i, let v f = (v\, vj) be the tangential and normal components of v e . What happens subsequent 
to impact depends on whether x e G dT Q or dD. In the case of a collision with 8T , the tracer 
bounces off dT with angle of reflection equal to angle of incidence, i.e., 

(v D e y = -vi (v\y = v\, (4.i) 

and the other variables are unchanged. In the case of a collision with the disk, the following energy 
exchange takes place between the disk and the tracer: 

(v D e )' = -vl (v\)' = u, J = v\. (4.2) 

Here we have, for simplicity, taken the radius of the disk, the moment of inertia of the disk, and the 
mass of tracer in such a way that the coefficients in Eq. (4.2) are equal to 1. The identification in the 
definition of fl k is z ~ z' where z,z' G M k are such that all of their coordinates are equal except 
that Vg and lo in z are replaced by the corresponding quantities with primes in Eqs. (4.1) and (4.2) 
for z' . We also write F(z) = z'. 

Note that in both (4.1) and (4.2), total energy is conserved, i.e., \v\ 2 + lo 2 = \v'\ 2 + (lo 1 ) 2 . The 
energy surfaces in this model are therefore 

n k:E = (r k xdDx sf)/ ~ 

where _ 

Sf = {(«!, . . . ,v k ,u) € R 2k+1 : \ v i\ 2 + " 2 = E}- 

We claim that the natural invariant measure, or Liouville measure, of the (discontinuous) flow 
<i> s on fl k is 

™k = (A 2 lr) fc X i"l\dD) X A 2fc+1 

where X d is d-dimensional Lebesgue measure and v d is surface area on the relevant d-sphere. Once 
the invariance of fh k is checked, it will follow immediately that the induced measures m k E = 
(A 2 |r) fc x (^llao) x v 2k on Cl k E are $ s -invariant, as are all measures on Cl k of the form ip(E)m k E 
for some ip : [0, oo) — > [0, oo). 

The invariance of fh k is obvious away from collisions and at collisions with dT . Because 
collisions occur one at a time, it suffices to consider a single collision between a single tracer and 
the disk. The problem, therefore, is reduced to the following: Consider $ s on and let M 1D 
denote the part of the collision manifold involving D. To prove that fh l is preserved in a collision 
with D, if suffices to check that for A C M 1D and e > arbitrarily small, 

We leave this as a calculus exercise. 
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4.1.2 Coupling to neighbors and heat baths 

We now consider a chain of N identical copies of the dynamical system described in Sect. 4.1.1, 
and define couplings between nearest neighbors and between end cells and heat baths. 

Let 7 L and 7 R be two marked subsegments of dT of equal length; these segments will serve 
as openings to allow tracers to pass between cells. For now it is best to think of T as having a 
left-right symmetry, and to think of j L and 7 R as vertical and symmetrically placed (as in Fig. 2), 
although as we will see, these geometric details are not relevant for the derivation of mean energy 
and tracer profiles. We call the segments j L and 7 R in the ith cell 7^ and 7 R l) . For i = 1, . . . , N — 1, 
we identify 7^ with 7 L * +1 \ that is to say, we think of the domains of the ith cell and the (i + l)st 
cell as having a wall in common, namely 7^ and 7l , and remove this wall, so that tracers that 
would have collided with it simply continue in a straight line into the adjacent cell. (See Fig. 2.) 




Figure 2: A row of diamond-shaped boxes with small lateral holes (made by removing vertical walls 
corresponding to 7 R = 7l +1) ) to allow the tracers to go from one box to the next. The shapes of 
the "boxes" can be quite general for much of our theory. The configuration shown is the one used 
in the simulations discussed in Sect. 4.5, but with larger holes for better visibility. 

Tracers are injected into the system as follows. Consider, for example, the bath on the left. We 
say the injection rate is q if at the ring of an exponential clock of rate q, a single tracer enters cell 1 
via 7^ . (Note that the rate q is not the injection rate per unit length of the opening 7^ but per unit 
time.) The points of entry and velocities of entering tracers are iid, the law being the one governing 
the collisions of tracers with 7^ That is to say, the point of entry is uniformly distributed on 7 L \ 
and the velocity v has density 

c e"^ 1 lull sin(y>)| dv , c= -£-=- , (4.3) 

where v G M 2 points into 7^ and 92 G (0, n) is the angle v makes with 7^ at the point of entry. 
(This is the distribution of v at collisions for particles with velocity distribution ^ exp(— (5\v\ 2 )dv.) 
Here (5 = 1/T where T is said to be the temperature of the bath. Observe that the mean energy of 
the tracers injected into the system by a bath at temperature T is not T but 3T/2. Injection from 
the right is done similarly, via the opening 7 R . When a tracer in the chain reaches 7^ or 7 R , it 
vanishes into the baths. 



Hamiltonian Models 



20 



This completes the description of our models. We remark that the process above is a Markov 
process in which the only randomness comes from the action of the baths. Once a tracer is in the 
system, its motion is governed by rules that are entirely deterministic. 

4.2 Single-cell analysis 

In analogy with Sect. 3.2, we investigate in this subsection the invariant measure for a single cell 
coupled to two heat baths with parameters T and g. 

Let Cl k and 0, kE be as in Sect. 4.1.1. As before, a state of this system is represented by a point 

in 

^ = U feLo^fc = U fc^0 U £>0 ^k,E 

where VL k and Cl k E are quotients of Q k and Q k E respectively obtained by identifying permutations 
of the k tracers. With { • • • } representing unordered sets as before, points in Q are denoted by z = 
({x 1 , . . . , x k }, {v x , . . . , v k },u>) or simply ({x e }, {v e }, uj), with v e understood to be attached to 
x t The quotient measures of fh k and fh kE are respectively m k and m kE . 

Abusing notation slightly, we continue to use $ s to denote the semi-flow on Cl, and let $ s 
denote the induced semi-flow on Q. Then <£ s is as in Sect. 4.1.1 except where tracers exit or enter 
the system. More precisely, if & s (z) G fl k for all < s < s , and a tracer exits the system at time 
s , then <E> S) (z) jumps to f^-i- Similarly, if a tracer is injected from one of the baths at time s , then 
instantaneously $ s (z) jumps to Q k+ i, the destination being given by a probability distribution. 

Let I7I denote the length of the segment 7 L or 7 R . 

Proposition 4.1 There is an invariant probability measure fi with the following properties: 
(a) the number of tracers present is a Poisson random variable with mean x where 

A 2 (T) q 



2^- 



(b) the conditional density of p, on Q k is c k a k dm k where 

<r k ({x t },0;{v e },u) = e-^+Stil-/) 
and c k is the normalizing constant. 

We observe as before that the Poisson parameter x is proportional to g (the higher the injection 
rate, the more tracers in the cell) and inversely proportional to y/T, i.e., the speed of the tracers (the 
faster the tracers, the sooner they leave). Unlike the models considered in Sect. 3, where the tracers 
are assumed to leave the cell at a rate equal to their speed, here the ratio A 2 (r)/|7| appears, as it 
should: the smaller the passage way, the longer it takes for the tracers to leave. 

Notice that we have not claimed that ^ is unique. 

We introduce some notation in preparation for the proof. For A C Q k and h > 0, we let 
&_ h (A) denote the set of all initial states in Q that in time h evolve into A assuming no new 
tracers are injected into the system between times and h. 1 Then &_ h (A) = U n>0 < ^ ) ^(^4) where 

Q_ h {A) = &__ h (A) n £l k+n , i.e. , Q^liA) is the set of states where initially k + n tracers are present, 
and by the end of time h exactly n of these tracers have exited and the remaining k are described by 
a state in A. 



'Notice that (1) h > 0} is a semi-flow, and is not defined; (2) <fr_ h (A) as defined is / (<E> h ) 1 (A). 
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Lemma 4.2 Let \i be as in Proposition 4.1, and let A £ be a cube of sides e intt k , e small enough 
that fi(A E ) « p k c k a k (z)e ik+2 for some z. We assume the following holds for all small h > 0: 

(i) no tracers are injected into the system on the time interval (0,h] ; 

(ii) <& h ($ ( °j! (A e )) = A £ . 
Then 

K*- h (A e )) = o k (z)e^ 2 [p k c k + 2h^ Pk+1 c k+1 + o(h)J (4.4) 
where p k = fJ<(£l k ) and c = -^L/3 3//2 . 

Proof: The idea is that for a particle to exit in the very short time h, it must be close to the exit 7 L or 
7 R and move towards it without colliding with the disk or the boundary, or it must have very large 
speed (and that is improbable). 

By assumption (i), we have fi(®_ h (A £ )) = Y, n >o M(^_^(^ £ ))- The n = term is handled 
easily: By virtue of (i) and (ii), the situation is equivalent to that in Sect. 4.1.1. Since fi\ n ^ is 

invariant for the closed dynamical system with k tracers, we have fj,($_ h (A £ )) = fJ,(A £ ). 

Consider next n = 1. We give the estimate for pL($P^ h (A e )) assuming 7 L and 7 R are straight-line 
segments, leaving the general case (where these segments may be curved) to the reader. First some 
notation: For v € M 2 and a > 0, let E(v , a,L) be the parallelogram on the same side of 7 L as V 
and with the property that one of its sides is 7 L while the other is parallel to v and has length a; 
E(v,a,R) is defined similarly. To simplify the discussion, we assume that for a > sufficiently 
small, E(v,a,L) and E(v,a,R) are contained in T, and leave to the reader the verification that 
"corners" at the end of 7 L or 7 R lead to higher order terms (in the variable h used below). 

Starting from a state in ^ } h (A £ ), we let x and v denote the initial position and velocity of the 
tracer that exits before time h, and treat separately the cases (1) \v\ < % and (2) \v\ > In Case 
(1), in order for the tracer to exit before time h, we must have x G E(v, h\v\ , L) U E(v, h\v\, R), and 
v must point toward the exits. Since X 2 (E(v , h\v |, L)) = X 2 (E(v, h\v\,R)) = h\v\ \ sin(</?)| \j\ where 
if is the angle v makes with 7 L or 7 R , we obtain 



-f3\v\ 2 



K^{(A £ )n{\v\ <^}) = a k (z)e 4k+2 p k+1 c k+1 -2h\ 7 \ T d^\ & m^)\ [ dv\v\e 
h Jo ^M<f 

= a k (z)e Ak+2 p k+1 c k+1 ■ 2%|i(l + o(h)) 

with c = 2(3 3/2 /^. For Case (2), we have the trivial estimate a k {z)e 4k+2 p k+1 c k+1 ■ o(h). 
To see that the terms corresponding to n > 1 are negligible, we first derive the bound 

< a k (z)e ik+2 Pk+n c k+n (2%|± + o(h)) n . (4.5) 
Then we compute the growth rate of p k+n c k+n . By the definitions of these numbers, we have 

\k + l Jm* J \k + l V | 7 | VTj 



Ck Pk 
giving 



>>!<+» '''■■+« " I ^ ) l} l-- r l-- " - ■ (4 - 6) 
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From (4.5) and (4.6) it follows that < a k (z)e ik+2 (const. ■ Kf. 

The asserted bound (4.4) for /j,(&_ h (A)) is proved. □ 

The main difference between the proofs of Propositions 3.5 and 4.1 is that Hamiltonian models 
have both geometry and memory. In preparation for the proof, we introduce the following language. 
Let A £ be as in Lemma 4.2. For £ = 1, . . . , k, we let denote the projection of A £ onto the plane 
of its a^-coordinate, and V £ the projection of A £ onto the plane of its ^-coordinate (so that X e and 
V e are e-squares in V and IR 2 respectively). We assume for simplicity that for each £, either X e is 
a strictly positive distance from 7 L and 7 R , in which case we say X e is in the interior, or one of its 
sides is contained in 7 L or 7 R . In the latter case, we say X e is adjacent to an exit. We further assume 
that if X e is adjacent to an exit, then either all v e <E V £ point toward the exit or away from it. 

Proof of Proposition 4.1: The invariance of \i is already noted in Sect. 4.1.1 except where it per- 
tains to entrances and exits of tracers. We focus therefore on these events, noting that the probability 
of more than one tracer entering on the time interval (0, h) is o(h), as is the probability of a tracer 
entering and leaving (immediately) on this time interval. These scenarios will be ignored. 

Let A £ be as above. We seek to show as before that J I A (z')P h (dz l \z)/j,(dz)\ h=0 = 0. Here 
it is necessary to treat separately the following configurations for A £ : 

Case 1. The following holds for all £: X^ can be in the interior or adjacent to an exit, and if it is 
adjacent to an exit, then all v £ in V t must point toward the exit. Notice that this configuration is 
relatively inaccessible, meaning the probability of a new tracer entering on (0, h) leading to a state 
in A £ is o(h)(i(A £ ). Notice also that this configuration has the property $ h (Q_ h (A E )) = A £ , so that 
the contribution of the no-new-tracers event to J I A (z')P h (dz'\z)p{dz) is, by Lemma 4.2, 

(1 - hgf a k (z)e 4k+2 ( Pk c k + 2h\j\^p k+1 c k+1 + o(h)) 

= a k (z)e ik+2 (p k c k (l - 2hg) + 2h\j\ ^P k+1 c k+1 + o{h)^j ( 4 . 7 ) 
= a k (z)e 4k+2 (p k c k + o(h)) , 



the last equality being valid on account of Eq. (4.6). 

Case 2. X 1 is adjacent to an exit and v 1 points away from it; X e and V e for £ > 1 are as in Case 
1. In this configuration, there is a part of X 1 that can only be reached in time h if one starts from 
outside. This region is a parallelogram similar to that in the proof of Lemma 4.2 but with one of its 
sides equal to X 1 Hj L or X 1 Dj R . Following the estimates in Case 1, we obtain that the contribution 
of the no-new-tracers event to J I A (z')P h (dz' \z)fi{dz) in this case is 

a k {z)e 4k+2 p k c k (l - sin(^)| + o(h)j (4.8) 

where v 1 is the v 1 coordinate of z and (p l is the angle v 1 makes with 7 L (or 7 R ). 

We now argue that the negative term above is balanced by the contribution of the event in which 
a new tracer enters on the time interval (0, h). This new tracer must have v 1 G V 1 and must enter 
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through the e-segment X 1 n 7 L or X 1 Pi j R . We claim that the probability of this event is 

Pk-iCfc-iafcC^e^ife 4 *- 2 • ^ • c| sin^lNe - ^" 1 '^ 2 . (4-9) 

The first factor in (4.9) is the ^-measure of the states corresponding to those in A £ but without the 
tracer with position and velocity (x ± , vf); the second factor is the probability of a tracer entering 
through the designated segment, and the third is the fraction of tracers entering with velocity G V l 
(see (4.3)). That (4.8) and (4.9) add up to n(A £ )(l + o(h)) again follows from (4.6). 

Case 3. X 1 and X 2 are adjacent to exits, v 1 and v 2 point away from the exits in question, and X e 
and V t are as in Case 1 for i > 2. We assume for simplicity that either (X 1 x V-f) n (X 2 x V 2 ) = 
or X 1 x V 1 = X 2 x V 2 . 

In the case (X 1 x V^) n (X 2 x V 2 ) = 0, the contribution of the no-new-tracers event is 

a k {z)e Ak+2 Pk c k (l - ^Nlsin^)! - ^ 2 ||sin(^ 2 )| + o(h)j , (4.10) 

and this is cancelled perfectly by the estimate corresponding to (4.9). 

In the case X 1 x V 1 = X 2 x V 2 , on Q k , where tracer positions and velocities are regarded as 
ordered /c-tuples, the set of states where both (x 1 ,v 1 ) and (x 2 ,v 2 ) are not reachable in time h is 
o(h), and the set where exactly one of these is not reachable is the union of two sets that project to 
the same set under ix k . Thus the estimates for both cases are as in Case 2. 

The remaining cases are handled similarly. □ 

Proposition 4.3 For the N -chain defined in Sect. 4.1.2 with T L = T R = T and g L = g R = g, the 
N-fold product fj, X • • • X fx is invariant. 

It suffices to check that the transfer of energy from one cell to the next leads to the correct 
relation between p k c k and Vk+\ c k+\- T ne proof is left to the reader. 

4.3 Derivation of equations of macroscopic profiles 

Having found the candidate family of Gibbs measures {fi T ' Q }, we now proceed as in Sect. 3.3, 
seeking to derive the relevant macroscopic profiles under Assumptions 1 and 2; see Sect. 3.3. There 
are two new problems, leading to two additional assumptions which we now discuss. 

The first problem is that of uniqueness and ergodicity. Unlike their stochastic counterparts, the 
Hamiltonian chains defined in Sect. 4.1 may not be ergodic; they are, in fact, easily shown to be 
nonergodic for certain choices of T . Without ergodicity, it is not clear how to make sense of the 
notion of local temperature, which lies at the heart of Assumption 2. Postponing a discussion till 
later, we bypass this issue by introducing 

Assumption 1'. We assume p, N is the unique invariant probability measure for the 
N -chain defined in Sect. 4.1. It follows that fi N is ergodic. 

Our next assumption pertains to the asymmetry of exit distributions from each cell in the finite 
chain. Let j N i denote the mean number of exits out of the ith cell per unit time in the iV-chain. 
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(Each time a tracer exits the 2th cell, it is counted as "an exit", even if it is the same tracer that 
re-enters and exits again.) Then j N i = j N i L + jjv,i,R where jjv j L anc ^ 3n,i,r 316 tne numbers of 
exits per unit time that go to the (i — l)st and (i + l)st cells respectively. In the random-halves model 
studied in Section 3, our local rule is that j Ni L = jw,i,R- That is typically not the case in models 
in which local rules are purely dynamical, such as the Hamiltonian models under consideration. 
Assumption 3 controls this deviation from symmetry for j N i and Q N i where j Ni is as above and 
Qn,i denotes the mean total energy transported out of the ith cell per unit time in the iV-chain. 

Assumption 3. For each cell configuration (TqjTljTr) and parameters T L , T R , q l , 
q r > 0, we assume 

(i) there exists an a > such that for all large N, the following hold for all i: 

1 . a 1 a 

\JN,i,R ~ 2~JN,i\ ^J-f i \Qn,i,R ~ gW^il < i 

1 1 a 

IO'jV,i,R ~ 2~jN,i) - (jN,i+l,R ~ 2~lN,i+d\ ^ ^2 ! 

1 la 

\(Qn,i,R - 2~QN,i) - (Qn,i+1,R - 2^N,i+l)\ ^ jyf ; 

(ii) as N —> oo, the profiles j N i and Q^,i tend to C 2 functions j and Q on [0, 1]. 



Remark 4.4 Assumption 3 is consistent with the following observations: For a cell in the iV-chain, 
the temperature difference between the cell on its left and the one on its right is of order \T L —T R \/N, 
so one expects the marginal of p, N at this site to deviate from the equilibrium measure in Sect 4.2 by 
the same order of magnitude. This deviation is in turn reflected in the differences \j N s jR — jj^ i L \ 
and \Q NiR — Qn,i,l\- Similarly, if the second differences are well behaved as we assume, their 
orders of magnitude as indicated above are dimensionally correct. Detailed dependencies of this 
asymmetry on the physical parameters are beyond the scope of this paper. 2 

We are primarily interested in situations where a ^ 1, which happens when there is good 
mixing within individual cells. Good geometry of T (such as concave walls and the absence of 
"traps") and small passageways between adjacent cells (so most tracers stay in the cell for a long 
time) are conducive to fast mixing within individual cells. The presence of large numbers of tracers 
also enhances mixing. 3 

Theorem 4.5 Under Assumptions 1, 1', 2 and 3, we have the following limiting profiles for the 
models in Sect. 4.1 as a in Assumption 3 tends to 0. 

• mean number of exits out of a site per unit time : 

M) = 2 (q l + (q r - qOO ; 

2 We thank H. Spohn for interesting correspondence on this point. 

3 It is important to distinguish between the following two levels of mixing: mixing within cells, and mixing in the 
chain. For example, small passageways between cells enhance mixing of the first kind but are obstructions to the latter. 
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• mean total energy transported out of a site per unit time : 

Q(0 = 3(&T L + ( QR T R - Qh T L )0 ; 

• mean stored energy at a site : 

(n _ i Q(Q i ar L + ( gR T R - 6l t l x 

U 3 j(0 2 g L + (qr ~ Ql)Z 

in the case g L = g R , this simplifies to s(£) = ^(T L + (T R — T L )£) . 

• mean tracer energy : e(£) = 2s(£) ; 

• mean number of tracers : 



A 2 (n pr- 

where I7I = |7 L | = |7 R | is the size of the passage between adjacent cells ; 
• mean total-cell energy : 

E(0 = s(0 + MOe(0 = s(0 + ^-V^(OM)- 

Proof: The proof follows that of Theorem 3.9. First we read off the pertinent information from 
Proposition 4. 1 for a single cell connected to two heat baths with parameters T and g: 

(i) stored energy has density -^2=e _/3y and mean s = \ ; 

(ii) tracer energy has density (3e~P x and mean T ; 4 

(iii) mean number of tracers, x = V^^y j 

(iv) mean total-cell energy, E = T(x + ^) ; 

(v) mean number of jumps out of cell per unit time, j = 2g ; 

(vi) mean total energy transported out of cell per unit time, Q = ^- • j = 3Tg . 

To prove (i), for example, we condition on the event that exactly k tracers are present. Integrating 
out all other variables, we obtain that the distribution of u is const. e 1 . Thus the distribution of 
s = uj 2 is as claimed. Items (ii) - (iv) are proved similarly, and (v) and (vi) are deduced from the 
fact that the cell is in equilibrium with the two baths. 

As in the proof of Theorem 3.9, we deduce next the profiles of j and Q. Here, 

JN,i,R ~ JN,i+l,L = JN,i-l,R ~ JN,i,L ■ 

This together with Assumption 3 implies that as N — > 00, the limiting profile j(£) has the property 
= 0{a). Since j(0) = 2g L and = 2g R , we have j(0 -► 2(g L + (g R - g L )0 as a -+ 0. 
Analogous reasoning gives Q(£) — > 3(g L T L + (g R T R — g L T L )Q as a — > 0. 

To deduce the remaining profiles, consider ^ G (0, 1). We let fi N — > p T( -& ,g ^\ and use the 
single-cell information above combined with the profiles of j and Q to identify T(£) and g(£). For 
example, the formula for s is obtained as follows: T = is from (v) and (vi), and s = \T is 

from (i). □ 



4 Note that this is the energy density when the tracers are in the box, to be distinguished from (vi). 
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4.3.1 Comparisons of models 

1. Predicted profiles for Hamiltonian and stochastic models. We observe that as functions, the 
predicted formulas in Theorem 4.5 are of the same type as their counterparts in Theorem 3.9 but 
the constants are different. The similarity stems from the fact that they are derived from the same 
general principles. The differences in constants reflect the differences in fi T,s , which in turn reflect 
the differences in local rules (see below). 

2. Relation between s and e. To highlight the role of the local rules in the profiles studied in this 
paper, we recall the relation between stored energy s and individual tracer energy in the various 
models encountered: 

(a) Random halves (Sects. 3.1-3.3): s = 2e. (At collisions, energy is split evenly on average, but 
the expected time for the next clock is longer for slower tracers.) 

(b) Stochastic models simulating Hamiltonian systems where both disk and tracer have a single 
degree of freedom (Sect. 3.6): s = e. 

(c) Hamiltonian models in which the disk has one degree of freedom and tracers have two (Sec- 
tion 4.1): e = 2s. 

To this list, we now add one more example, namely 

(d) Hamiltonian models in which the disk has one degree of freedom and tracers have 3: Consider 
the model described in Sect. 4.1, but with T C M 3 and the disk replaced by a cylinder 
that rotates along a fixed axis. Here, Liouville measure for a closed system with k tracers is 
fh k = (A 3 | r ) fc x {v\\qd) x A 3fe+1 (cf. Sect. 4.1.1). From a single-cell analysis similar to that in 
Sect. 4.2, /j T ' e is easily computed. One notes in particular that the distribution of tracer energy 
is const. \fxe~^ x , while disk energy is as before. A simple computation then gives e = 3s. 

These examples demonstrate clearly that the relation between s and e is entirely a function of the 
local structure. In the case of Hamiltonian systems, we see that it is also dimension-dependent. 

The interpretation of results in Sect. 3.5 applies to all of the models above. 
4.4 Ergodicity issues 

Questions of ergodicity for the chains in Theorem 4.5 are beyond the scope of this paper. We include 
only brief discussions of the following three aspects of the problem: 

1. Randomness in the injection process 

Among the various features of our models, the one the most responsible for promoting ergodicity 
is the randomness with which new tracers are injected into the system. We observe, however, that 
this genuinely stochastic behavior occurs only at the two ends of the chain, and even there, the 
transition probabilities do not have densities with respect to the underlying Lebesgue measure. The 
problem is thus one of controllability involving the deterministic part of the dynamics. 

2. Hyperbolicity of billiard dynamics: a necessary condition 

Let A N C M? denote the playground for the tracers in the iV-chain. That is to say, it is the union 
of N copies of F arranged in the configuration shown in Fig. 2 with open passages between adjacent 
copies of T. The presence of one of more tracers being trapped in A N without contact with any 
of the turning disks or the openings at the two ends (i.e., 7^ and 7^) is clearly an obstruction to 



Hamiltonian Models 



27 



ergodicity. This scenario is easily ruled out by choosing T to have concave (or scattering) walls. 
Such a choice of T implies that dA N also has concave boundaries, and the free motion of a particle 
in a domain with concave boundaries is well known to be hyperbolic and ergodic [18, 13]. 
We do not know if the absence of trapped tracers in the sense above implies ergodicity. 

3. Enhancing ergodicity 

Without (formally) guaranteeing ergodicity, various measures can be taken to "enhance" it, 
meaning to make the system appear for practical purposes as close to being ergodic as one wishes. 
For example, one can introduce more scattering within each cell by increasing the curvature of the 
walls of T , or alternately, one could add convex bodies inside T that play the role of Lorentz scat- 
terers. Another possibility is to add a small amount of noise, and a third is to increase the injection 
rates: physical intuition says that the larger the number of tracers in the system, the more likely 
stored energy will behave ergodically. 

4.5 Results of simulations 

To check the applicability of the theory proposed in Sects. 4.1-4.3 to real and finite systems, we 
have done extensive simulations some of which we describe in this subsection. The domain T used 
in our simulations is as shown in Fig. 2. Actual specifications of T are as follows: We start with 
a square of sides 2, subtracting from it first 4 disks of radius 1.15 centered at the 4 corners of the 
square. Two openings corresponding to 7 L and 7 R are then created on the left and right; each has 
length 0.02. This completes the definition of T . The disk D is located at the center of the square; 
it has radius r = 0.0793. 

Our choice of domain was influenced by the following factors: First, dT is taken to be piece- 
wise concave to promote ergodicity. Second is the size of the disk: A disk that is too small is hit by 
a tracer only rarely; many tracers may pass through the cell without interacting with the disk (this is 
analogous to having a large 5 in Sect. 3.1). A disk that is too large (relative to the domain in which 
it can fit) may cause an unduly large fraction of tracers entering the box to exit immediately from 
the same side. Both scenarios lead to large time-correlations, which are well known to impede the 
speed of convergence to n N in a finite chain. They may also affect the infinite- volume limit. 

We have found the geometry and specifications above to work quite well, with a tracer making, 
on average, about 71 collisions while in a cell. Of these collisions, about 12.5 are with the disk. 

For the single cell (with the geometry above) plugged to two identical heat baths, we have tested 
the system extensively for ergodicity. To the degree that one can ascertain from simulations, there is 
an ergodic component covering nearly 100% of the phase space. The various energy distributions as 
well as the Poisson distribution of the number of tracers present agree perfectly with those predicted 
by Proposition 4.1. 

Simulations for chains of 20 to 60 cells with the choice of r and |— y| above showed very good 
agreement with the theory. A sample of the fits for Q(£), s(£) and E{£) for 9 • 10 9 events and 30 sites 
is shown in Figure 3. Here the bounds in Assumption 3 are very small, that is to say, the ejection 
rates to the left and right are very close to 50/50. We have also investigated the quantity a in 
Assumption 3 for various values of r and |— y|, up to r = 0.23 (which is quite close to the maximum- 
size disk that can be fitted into the domain T ) and I7I = 0.06. Our findings are consistent with the 
discussion in Sect. 4.3. 

In addition to these profiles, we have also verified directly Assumption 2, which asserts that the 
distributions of energy and tracers within each cell are in accordance with those given by \F ,e for 
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some T, g depending on the cell. A sample of these results is shown in Fig. 4. 
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Figure 3: Rotating disks model with chain of 30 cells, temperatures T L = 100, T R = 10, and 
injection rates g L = 1, g R = 2. Top left: Q t , energy transported out of site i per unit time as a 
function of i. Top right: Mean disk energy Sj. Bottom: Mean total energy 



4.6 Related models 

In this subsection we recall from the literature a few models that in their original or slightly modified 
form can be regarded as approximate realizations of the class described in Sect. 2. 1 of this paper. 
For more complete accounts, see the review papers [1, 10, 12]. 

The models which come closest to ours, and which to some degree inspired this work, are those 
in [14, 9]. In these papers, the authors carried out a numerical study of a system comprised of an 
array of disks similar to those in Sect. 4. 1 but arranged in two rows with periodic boundaries (in the 
vertical direction). These disks interact via tracers following the rules first used in [17]. We have 
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Figure 4: Same parameters as in Fig. 3. Top 2 figures show semi-log plots of tracer energy distribu- 
tions at various sites. Top left: Densities of tracer energies inside boxes (theory predicts f3e~@ x ). Top 
right: Densities of tracer energies upon exiting the various boxes (theory predicts 2/32 ^x/ire~P x ). 
Bottom: Distribution of numbers of tracers at several sites (theory predicts Poisson distribution). 
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adopted the same local rules, but have elected to arrange our disks in a single row to simplify the 
analysis. 

There is a number of papers dealing with mechanical gadgets that on some level appear similar to 
ours. For example, in [1 1, 7], vertical plates are pushed back and forth by particles trapped between 
them. The main difference between these models and ours is that they have exactly one "tracer" 
in each "cell". In this respect, these models are closer to our earlier work [5] in which locked-in 
tracers were considered. Ding-a-ling and ding-dong models belong essentially to the same class 
[2,16,6,15]. 

We mention that nonlinearities of profile are difficult to see when the temperature differences at 
the two ends are relatively small (in fact, what counts in many cases, including the models studied 
in this paper, is the ratio of temperatures at the two ends). This may explain why some authors have 
reported linear profiles when our analysis suggests that may not be the case. 

Finally, we mention a very well-studied situation, namely that of the Fermi-Pasta-Ulam chain. 
In this model, and in many others, there is a potential of the form 

V(x i - x i+1 ) + ViXj) , 

with U and V functions that grow to oo and x i the coordinates of a chain of anharmonic oscillators. 
The pinning potential V plays the role of the "tank" in our models, while the interparticle potential 
is more akin to the role of the tracers. This class of models is difficult to handle because in contrast 
to the basic setup in our study, there is no clear separation of the pinning and interaction energies. 

Acknowledgments: The authors thank O. Lanford for helpful discussions. JPE acknowledges the 
Courant Institute, and LSY the University of Geneva, for their hospitality. 
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